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NOMENCLATURE 


SYMBOLS 

B(  )  correlation  function 

(  )  correlation  function 
second  derivative 

c  acoustic  phase  velocity; 

sound  speed 

E  expected  value 

Af  doppler  frequency 

(Appendix  II) 

J  Jacobian  (Eq.  30) 

N  number  of  multipaths  or 

arrivals 

p(*)  probability 

r  correlation  distance 

X  distance  downrange  from 

source 

y  distance  crossrange, 

normal  to  x-z 

2  distance  above  reflecting 

bottom 

0  bottom  slope 

CT  standard  deviation 

V  speed 

cp  acoustic  ray  angle  to 

horizontal 

At  pulse  time 


SUBSCRIPTS 

I  lower 

M  multipath 

max  maximum 

min  minimum 

o  surface  to  bottom  or  reference 

value 

R  receiver/reflector 
S  source 

R  round-trip 
u  upper 

z  elevation 

0  slope 

1  outgoing  or  one-way 

2  returning 

c  center 
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I.  INl’FDDUCTION 


The  reflection  of  acoustic  energy  by  a  rough  suj  Pace,  such  as  the  ocean 
bottom,  is  of  practical  interest.  If  the  acous  Ic  signals  propagated  in 
such  a  manner  occur  between  moving  source  and  receiver  or  reflector, 
and  if  the  reflecting  boundary  is  not -too -rough,*  then,  due  to  a  single 
pulse,  there  may  be  numerous  more  or  less  discrete  "arrivals”  (Refs, 

1,2)  at  both  terminals  of  the  propagation  path.  These  arrivals  which 
may  have  an  associated  doppler  shift  in  frequency  arrive  along  various 
"multipaths."  In  signal-processing  techniques  using  time  and  frequency 
correlation,  an  understanding  of  relationships  among  kinematic,  geo¬ 
metric,  and  environmental  aspects  of  such  a  problem  is  essential.  This 
paper  presents  an  analysis  for  detexirining  on  a  statistical  basis  ex¬ 
pected  acoustic  sigial  arrival  times  and  corresponding  expected  doppler 
frequency  shifts  for  both  one-way  and  round-trip  paths.  There  may  also 
be  an  ionospheric  radar  application. 

This  paper  proceeds  by  analyzing  the  transit  time  and  doppler 
frequency  characteristics  of  a  narrow  pulse  on  one-way  multipaths  from 
a  sea  bottom  with  random  normal  distribution  of  slope  but  negligibly 

V/  ■ 

A  not -too- rough  surface  is  one  which  has  variances  of  surface  eleva¬ 
tion  derivatives  sufficiently  small  that  there  is  no  appreciable 
occultatlon  of  surface  area  at  oblique  incidence  angles;  as  a  related 
practical  matter  in  transmitting  energy  to  a  receiver  or  reflector, 
it  is  necessary  that  surface  correlation  lengths  be  sufficiently 
great  that  at  least  the  first  Fresnel  zone  is  nominally  flat. 
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large  bottom  elevation  variance;  it  continues- -through  a  convolution 
of  this  "one-way”  case — by  determining  the  properties  of  narrow  pulses 
on  round-trip  multipaths;  next,  this  round-trip,  narrow  pulse  analysis 
is  extended  to  long  pulses  anu,  finally,  the  importance  of  non- 
negligible  bottom  elevation  variance  and  source  and  receiver-reflector 
motion  is  considered. 


II.  ANALYSIS 


A.  ONE-WAY  PATHS,  MRPDW  PULSE 

Let  Fig.  1  represent  one  of  a  number  of  one-way  paths  between  a 
moving  acoustic  source,  S,  and  a  moving  receiver,  R,  i.e.,  a  narrow 
pulse  transmitted  from  S  is  r'":ived  at  R  inteii<-lunally  (communicating) 
or  unintentionally  (interferring) .  The  analysis  is  restricted  for  now 
to  narrow  (delta-function)  pulses  for  simplicity  but  the  more  general 
finite  width  pulss  will  be  treated  subsequently  with  a  measure  of 
"narrow"  given  then.  In  this,  0  is  the  local  slope  of  the  reflecting 
boundary  which  Is  presumed  to  have  correlation  lengths  of  bottom  ele¬ 
vation  and  its  first  two  derivatives  very  large  relative  to  the  wave¬ 
length  of  the  acoustic  energy.  As  a  concomitant  of  this,  it  is  pre¬ 
sumed  also  that  bottom  elevation  variance  is  unimportant.  As  this 
latter  assumption  vd.ll  not  always  be  valid,  the  effect  of  importantly 
large  elevation  variance  1  ^  treated  in  a  subsequent  section.  For  th® 
sake  of  simplicity  in  disc  ssing  the  problem,  only  one -dimensional 
roughness  in  the  direction  of  propagation  will  be  considered;  extension 
to  two-dimensional  roughness  is  possible  though  tedious  (Appendix  I). 
The  one-dimensional  problem  is  reasonably  descriptive  of  the  tvio- 
dimensional  situation  because  the  probability  of  paths  at  large 
azimuthal  excursions  from  the  plane  of  Fig,  1  is  small  (Ref.  2)  for 
typical  sea  bottom,  and  this  being  so,  the  signal  transit  time  and 
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dopplar  frequency  between  S  and  R  are  determined  with  only  small  error 
by  in-plane  components. 

Now  in  Fig.  1,  the  relationships  oL.iong  the  several  variables  is 
as  follows:* 


X  =  Zs  coc  cp 


(lA) 


x*  -  X  =  Z(,  cot  (cp  -  26) 


(IB) 


from  which  it  develops  that  a  ray  from  S  will  reflect  to  R  if 


M 


Cp  -  arccot 


Zs  cot  cp 

Zn 


)] 


(2) 


where  the  subscript  M  indicates  that  a  multipath  exists  between  S  and 
R  as  a  function  of  cp,  z, ,  Zs ,  and  x* .  It  has  been  shown  (Refs.  1,2; 
Appendix  I)  that  if  6  has  a  Gaussian,  l.e.,  normal  distribution  then 
the  probable  or  expected  number  of  multipaths  in  an  increment  of 
bottom  distance  x  is 


p(M  in  dx:  cp,  z*  ,  Zr  ,  x* )  »  dE(N^ )  =  - 


exp  (-6m/2oJ)  dx  (3) 


It  may  be  mathematically  more  satisfying  to  assume  tar  9  as  having 
Gaussian  distribution;  if  cTq  «  1  it  makes  little  difference  which 
assumption  is  made,  although  as  9  s  tan  0,  Gaussian  6  yields  greater 


•k 

If  index  of  refraction  variation  is,  on  the  average,  important  then 
this  analysis  must  be  modified  (Ref.  2). 
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expected  values  as  is  clear  from  Eq.  3.  Presently,  experimental  data 
are  insufficient  to  make  a  decision  between  the  two  assumptions  and  a 
Gaussian  oistribution  of  9  is  used  to  simplify  calculation. 

Nov;  the  6"'pected  niimber  of  multipaths  in  dx  may,  using  Eqs.  lA 
and  3,  be  transformed  to  the  expected  angular  density  of  one-way  multi¬ 
paths  E(N^  )  as 

[-  -  -P 

In  anticipation  of  determining  transit  time  and  doppler  frequency 
of  multipath  arrivals,  the  cumulative  expected  nurnher  of  multipaths 
E*(cp)  between  cp  =  0*  deg  and  (p  is  determined  from  Eq.  4  as 

[^] 

tha  integration  being  carried  from  <p  =  O'*’  to  avoid  a  singular  point  at 
cp  =  0  which  arises  from  using  0  rather  than  tan  0  in  Eq.  3. 

Given  the  conditions  described  graphically  in  Fig.  1,  signal 
transit  time  and  doppler  fi-equency  may  be  determined.  From  Fig.  1  the 
one-way  signal  transit  time  of  narrow  pulses  is 

t^/(zs/c)  =  ^csc  cp  +  (Zn/Zj)  CSC  (cp  -  29)]  (6) 
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and  the  one-way  doppler  frequency  is,  approximately  (Appendix  II) 


(Afi  /fo)~ 

Vs/c 


2:  COS  Cp  +  (Vr/Vj)  cos  (cp  -  26) 


(7) 


In  Eqs.  6  and  7,  c  is  sound  speed,  fo  the  at-rest  frequency  of  the 
sound  propagation  of  S  and  it  is  assumed  that  Vj  and  Vr  are  positive 
as  shown  in  Fig.  1.  Figure  2  shows  E*(cp),  t^  (cp)  and  Af (cp)  for  the 
following  values  of  para'neters  to  be  used  in  this  paper  unless  other¬ 
wise  noted: 


Z5  =  Zr 
Vj  =  Vr 

c 

% 


=  5.5 


2200  fathoms  =  13,200  ft  =  4024  m 
0.01  c 

5000  fps  =  2961  knots  =  1524  m/sec 
arccot  [xr/Czr  +  Zr)]  =  20  deg 
15  deg 
25  deg 

3.75  deg  -  0.0655  radian 

_  /-inoc  —  ^■*10  c 


(8) 


[-BJ(0)/Bq(0)]^  =  (1025  ft)"‘  =  (312.5  m)" 


Figure  2  has  special  utility  in  that  if  transmission  from  S  is  between 
cf^in  and  Cf^Rx  ^nly,  then  for  this  beamwidth  the  expected  number  of  one¬ 
way  multipaths  between  S  and  R  is  just  E*(C(\irx)  " 
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FIGURE  2  One-Way  Multipath  Expected  Number,  Transit  Time  and  Doppler 
Frequency— Narrow  Pulse 


1.  Multipath  Maximum  Number  for  Fixed  Beamwidth 

Equation  5  which  gives  an  expression  for  expected  number  of  multi¬ 
paths  may  be  used  to  determine  beam  depression  to  maximize  the  number 
of  arrivals  at  R  from  S  for  a  given  beamwidth,  Acp.  Consider 


(^)  ['wJ  A- 


•%+iAcp 
'cPe  -iACp 


CSC  cp  exp 


[-(9S/2c;)] 


dcp  (9) 


where  Acp  is  a  fixed  beamwidth  and  cpe  is  the  beam  center.  If  Eq.  9 
be  differentiated  with  respect  to  cpc  and  set  equal  to  zero,  then  there 
results  the  condition  for  which  E(Nj^ )  is  maximized  as  a  function  of 

Using  Leibnitz  rule,  the  derivative  of  Eq.  9  with  respect  to  % 

is 


dE(Ni  ) 

"3^ 


(lOA) 


CSC*  (cp,  -  lAcp)  exp  [^-eM(cpe  -  iAcp)/2a^j  +  o| 


If  this  derivative  be  set  equal  to  zero  there  results 


CSC®  (cpo  *■  iAcp)  exp  ^-0^(Cft,  + 


(lOB) 


CSC®  (cp,  -  ^Acp)  exp  ^-0S(cft,  -  iAcp)/2<3^J 
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which  by  reference  to  Eq.  4  is  simply  the  condition  that 


■^1  = 


.  In  other  words,  for  a  given  beamwidth  -  %iii  -  Acp, 

the  number  of  multipaths  within  Acp  will  maximize  when  the  angular  den¬ 
sity  of  multipaths  at  the  beam’s  edges  are  equal;  generally  cpc 
arccot  [X|,/(zs  +  Zr)]  =  %>  the  apparent  specular  point. 

2.  Doppler  Frequency  Mean-Slope  Deviation  Relationship* 

A  measure  of  surface  slope  standard  deviation  from  the  doppler 
frequency  mean  devolves  from  the  case  in  which  S  follows  R  at  fixed 
X|, ,  i.e.,  |vs|  =  |vi,|.  Then,  for  9  small,  doppler  frequency  is  given 
by 


cos  cp  -  c  '  (cp  -  20) 


—  -  20  sin 


(11) 


and  for  this  condition  the  expected  value  of  the  absolute  deviation  of 
doppler  shift  is,  for  Og  small  also 


29  sin  cp 
00 


exp  (-0*/2o^)  d0 


(12) 


It - 

What  follows  may  bear  some  pertinence  also  to  an  electromagnetic 
source  and  its  receiver  moving  relative  to  a  randomly  refracting 
ionosphere. 
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where  cp  =  cp(9)  from  Eq.  2.  If  a  mean  value  of  cp  s:  cpo  =  arccot 
x»/(zs  +  z, )  io  used,  then 


E 


(Af 1 /fo ) 

(vIM" 


40^  sin  % 


(13) 


i.e. , 


E|(Afi)|  s:  2  (vs/c)foaQ  sin  cpo  (14) 


Apparently,  therefore,  mean  doppler  shift  is  a  direct  measure  of  ocean 
bottom  roughness  variance;  for  conditions  of  Eq.  8,  Eq.  1.  takes  a 
value  of  0.2  cps/kcps.  A  related  result  is  obtained  if  <(Afi)®>  is 
determined  instead  of  < | Af x | ^  . 

3.  Transit  Time  and  Doppler  Frequency  Characteristics 

Figure  2  and  the  analysis  leading  to  it  provide  a  basis  for  de¬ 
scribing  some  of  the  statistical  characteristics  of  signal  transit 
time  and  doppler  frequency.  Generally,  the  expected  number  of  arrivals 
until  ti^  is 


where  cp,  (t/ )  >  cft,  and  cp^(tx' )  <  cft,  are  angles  appropriate,  as  in  Fig.  2, 
to  a  transit  time  t/ .  Thus  the  time  rate  of  arrivals  is  just  the 
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derivative  of  Eq.  15  which)  using  Leibnitz  rule,  is  (with  primes 
omitted)* 


dE(Ni  ) 
TtT^ 


!(^l  -li^l («■)!, 


cp#(t:  ) 


(16) 


where  cp^m  cp^jCti  )  <  cp,  and  cpo  <  cpu  (ti  )  ^  %iax  >  %  =  arccot 
x,/(zs  +  Zb);  if  cp^(ti)  <  or  if  cp^  (ti  )  >  cp,,,  then  dcp/dtj  -  0. 

The  expected  number  of  arrivals  from  S  at  R  prior  to  a  given 
time  t/  is  available  in  an  alternative  form  from  Eq.  16  by  integration; 
thus 


>]t.  *  t,' 


(17; 


where  the  integrand  of  Eq,  17  comes  from  Eq.  16  with  the  associated 
restrictions,  and  t,!#  is  the  least  transit  time  appropriate  to  the 
bearr,,  (cpkax,  %i,).  Clearly 


lim 

tj'-^ 


dE^ 

dti 


dtj 


(18) 


A  parallel  construction  for  doppler  frequency  characteristics  is 
possible  and  one  may  write  straightforwardly  the  analog  of  Eq.  16  as 


It  is  assumed  that  dX|,/dt  =  0  as  a  consequence  of  small  Vs  and  v* ; 
this  introduces  a  negligible  error. 
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(19) 


dE(Ni 


% (Afl  ) 


where  %in  ^  cp^CAfi  )  <  cp/  and  cp/  <  cp^CAfi )  ^  %»x  •  In  this  cp^  = 
tp[(Afi),xx]  from  Eq.  7;  again  if  cp^  and  cPu  lie  outside  the  interval 

(%  1  n  >  %  a  X  )  then  ^  ^  • 

The  expected  number  of  arrivals  from  S  at  R  with  doppler  frequency 
less  than  a  given  value  Afi^  is,  in  parallel  with  Eq.  17, 


[='"■>] if,  .if,' 


(20) 


The  transit  time  and  doppler  frequency  characteristics  of  the 
signals  are  not  independent  as  Fig.  2  suggests;  For  a  given  one-way 
path  arrival  time,  one  of  two  possible  sonar  beam  depression  angles  is 
permissible,  in  general,  and  the  doppler  frequency  is  fixed  by  these 
angles.  Hence  for  one-way  trans-mission,  doppler  frequency  is  one  of 
two  values  fixed  by  arrival  time. 

Figure  3  shows  the  cumulative  number  of  "arrivals,”  i.e.,  multi¬ 
paths  of  a  single  very  narrow  pulse  as  a  function  of  time  and  doppler 
frequency.  For  the  typical  values  of  Eq.  8,  numerical  values  from 
Fig.  2  are 

maximum  number  of  multipaths  or  "arrivals"; 

[E(Ni)j  =  1.32  ^13,200  ft^  ft^  ^  =  5.4 
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FIGURE  3  One-Way  Multipath  Transit  Time  and  Doppler  Frequency  Cumulative 
Density — Narrow  Poise 


and  from  Fig.  3, 


time  interval  of  50  percent  of  the  arrivals: 

(At)o.5  =  (5.854  -  5.848)  '^50§g°fps  '=  15.8  msec 

Figure  2  indicates  that 

maximum  doppler  frequency  (Afi ; 


18.8  cps/kcps  of  fo 


and  Fig  3  shows 

frequency  interval  of  90  percent  of  the  arrivals; 

A  =  (1-880  -  1.8535)  =  0.265  cps/kcps  of  f, 

i.e.,  substantially  all  the  arrivals  have  doppler  frequercy  between 
18.5  and  18.8  cps/kcps  for  R  closing  on  8  each  at  one  percent  sound 
speed,  and  for  typical  geometry  and  ocean  bottom  conditions. 

The  situation  which  Fig.  2  typifies  may  be  elaborated.  If  S  Is 
not  transmitting  omnidirectionally  (0  ^  <p  <  180  deg)  but  is  limited 
(<¥kin  then  the  expected  number  of  multipaths  to  R  may  as 

stated  be  determined  as  the  difference  EfCcfy,,)  -  E*(%ib)  for  this 
latter  case.  The  distributions  of  Fig.  3  will,  of  course,  be  modified 
by  this,  but  these  actual  distributions  may  be  determined  from  Fig.  2 
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with  care  for  beam  limits.  If  further,  R  is  "listening”  with  a  limited 
beam  then  the  number  of  multipaths  and  their  time  and  frequency  charac¬ 
teristics  will  be  determined  by  the  interval  Ax  of  ocean  bottom  which 
is  common  to  the  beamwidths  of  both  S  and  R.  Finally,  if  the  upper 
bounding  surface  is  planar  for  waves  of  frequency  fo  (Ref.  3)  then  the 
number  of  paths  to  R  from  S  via  the  surface  may  be  determined  by  re¬ 
placing  S),  in  the  foregoing  by  2zo  -  (the  distance  above  the  bottom 
of  the  virtual  receiver).  If  (Zo  -  z,,  )/Zo  «  1  then  the  number  of 
paths  via  the  bottom  and  surface  is  approximately  that  via  the  bottom 
alone,  and  the  arrival  time  and  doppler  frequency  characteristics  will 
be  similar  also  for  the  two  path  types.  Hence  for  this  situation, 
numbers  of  arrival  and  densities  of  arrivals  in  time  and  frequency  will 
be  approximately  doubled. 

B.  ROUND-TRIP  PATHS,  NARROW  PULSE 

The  treatment  above  of  one-way  paths  provides  the  basis  for  deter¬ 
mining  the  expected  nur.iber,  transit  time,  and  doppler  frequency  of 
round-trip  paths  for  sound  waves  reflected  from  a  not-too-rough 
boundary.  Figure  4  gives  the  geometry  and  nomenclature  for  this  round- 
trip  case  which  is  a  straightfor^;ard  elaboration  of  Fig,  1. 

Now  if  p(M  in  dxi  :  cpj  ,  2$ ,  z^ ,  x,)  is  the  probable  number  of  out¬ 
going  multipaths  occurring  in  dxi  ,  and  p(M  in  d)^  j  cffe  ,  z, ,  z» ,  •/, )  is 
the  corresponding  number  of  a  return  path  occurring  in  dxa  ,  then  the 
probable  number  of  round-trip  paths  as  a  result  of  reflection  at  R 
going  through  dxi  and  dxg  is  the  product  of  these  individual  incre¬ 
mental  niimbers  and 
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d^ECHa) 


;^]  [^]  -- 


(21) 


The  expected  number  E(l^  )  of  round-trip  paths  between  S  and  R  is  the 
integral  of  Eq,  21,  i.e,, 


E(Jt ) 


d^ECrt ) 


(22) 


or  in  explicit  form  as  a  function  of  beam  depression  angle 


E(I^) 


/,  \a  r  b:(o)’|  /-(cpi)  ■  ax  /•(4fe  )  ■  a  ] 

1,.. 


exp 


CSC  cpi  CSC  cfla 

(23) 


As  cpt  and  cc^  and  independent,  if  R  is  ensonlfied  and  reflects  only  in 
the  two  downward  quadrants  and  if  the  transmitting  and  receiving  beams 
at  S  are  identical,  then  from  Eqs.  4  and  23 


E(H»)  =  [E(Ni)p 


(24) 


If,  however,  R  is  ensonified  from  both  planar  surface  and  rough  bottom 
and  can  reflect  these  arrivals  in  either  the  bottom  or  surface-bottcxn 
paths,  then  the  tf'tal  expected  number  of  multipaths  may  be  as  large  as 
4rE(Ni)]®.  For  the  example  of  the  preceeding  section  for  one-way  paths 
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with  Zr  2s  Zo  ,  this  takes  a  value  of  4(5.4)®  >  116,  not  to  be  sure,  all 
of  the  same  intensity.  Figure  5  shows  the  density  of  multipaths  in 
transmit-angle  (cpj  )  and  receive-angle  (%  )  space  based  on  the  inte¬ 
grand  of  Eq.  23  and  it  is  interesting  that  the  maximum  angular  density 
occurs  for  cp  2:  18.7  degrees,  at  lesser  depression  than  the  "specular” 
angle  cp^  =  20  degrees.  Figure  5  is,  of  course,  a  two-dimensional 
representation  of  a  surface  and  provides  a  basis  for  considering  the 
transit  time  and  doppler  frequency  characteristics  of  round-trip 
multipaths  and  their  density  distribution  in  time- frequency  space. 

The  round-trip  transit  time  t,  of  a  narrow  pxase  is 

t(,  =  ti  +  ta  (25) 


where 


ti/(Z8/c)  =  [esc  cp,  +  (Zr/Zs)  CSC  (<p|  -  261)],  1  =  1,2 

Thus  once  the  geometry  of  Fig.  4  is  fixed,  ti  *  ti(<Pi).  The  round- 
trip  doppler  frequency  is 


Af,  =  Afx  +  Afa  (26) 

where  there  is  associated  with  each  direction  a  doppler  frequency 
given  approximately  as  (see  Appendix  II) 
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Receive  beam  ray  depression  angle,  (deg) 
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Tranunit  b«am  ray  d«praMton  angit,  (d«g) 


FIGURE  5  Round-Trip  Multipath  Dansity  in  Tranimit— and  Racaiva-Boam  Anglo 
Space— Narrow  Pulie 


Rt-n-M-S 

20 


(Afj^ 

(Vs/C) 


=  [^cos  cpi  +  (v,/vs)  cos  (cpi  -  2ei)J,  i  =  1,2 


and  again  =  Af i  (cpj ) . 

Given  the  foregoing  equations,  it  is  of  interest  at  first  to  de¬ 
scribe  the  arrival  rate  of  a  narrow  pulse  on  round-trip  paths.  Now 
the  probability  of  an  arrival  during  dt*  is  just  the  probability  of  an 
arrival  at  R  during  dt^  multiplied  by  the  probability  of  an  arrival  at 
S  from  R  during  dtg  =  d(t*  -  t^),  summed  over  all  t*  =  t^  +  tg  .  This 
is  shown  schematically  in  Fig.  6.  Thus 


p(M  in  dt* 


=  r  ■ 

•'(ti  ).i 


-  (to  )■  ll 


p(M  in  dtx )  p(M  in  dt«  )  dti  (27) 


or,  omitting  primes 


/dE(l^)\  _ 

i^tg  -  (ta)BlB 

’dE(Ni ) 

dEtKL} 

•'(t,).,. 

(H:, 

u  • 

t, 

Ttr^ 

m  V 

dt, 


Jtn  -  ti 


(28) 


with  the  constraint  that  all  of  the  Interval  ti ,  te  lies  within  the 
interval  •  Considering  the  complexity  of  the  relationships 

among  t^ ,  tg  and  cp^ ,  ,  It  Is  likely  that  the  maximum  of  Eq.  28  Is 

most  readily  found  through  numerical  integration  of  that  equation.  A 
construction  similar  to  Eqs.  27  and  28  yields  for  dopper  frequency 


‘dE(lfe  ) 

• 

^(Afi)..x 

dE(NiJI 

dECN, ) 

dAf, 

Afu  -  (Afa  )«a  X 

dAf, 

(Af,) 

dAlI 

dAf  (29) 
Af, 
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FIGURE  6  Schematic  for  Round-Trip  MulHpoth  Trontit-Timo  Gilculotion- 
Nkirrow  PuIm 


and  the  maximum  of  this  may  again  be  found  numerically.  Figure  7 
(obtained  by  graphical  calculation  and  subject  therefore  to  some 
error)  shows  the  time  and  frequency  distributions  resulting  from  Eqs. 
28  and  29.  For  the  parametric  values  of  Eq.  8,  the  interesting  points 
of  Fig.  7  are: 

transit  time  tn  for  maximum  arrival  rate:  (tn),>x  ran 

r*te  -  (t,  ),ia  ^  (11.699  -  11.695)  (13,200/5000)  =  16  msec 

maximum  rate  of  arrivals: 


s:  58  [(13,200)(5000)/tt®]  (1025  ft)"*  =  0.37/msec 

doppler  frequency  corresponding  to  m'*xiir.um  spectral  density: 

4^  as  3.757  (50  fps/5000  fps)  «  37.57  cps/kcps  of  fo 
maximum  spectral  density  of  arrivals; 

m  (It)  (^°°° 

The  time  and  frequency  distributions  of  Fig.  7  are  each  without 
relation  to  the  other.  Thus  one  may  desire  the  arrival  density 
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distributed  in  both  time  and  fi'equency.  Equation  21  may  be  trans¬ 
formed  from  (Xj ,  Xg  )- space  to  (t,, ,  Af(,)-space  through  the  Jacobian; 
thus 


d^E(^fc) 


J 


(Xi  ,  Xg  ) 
(t«,  Af,) 


dt(,  dAfj, 


where 


(30) 


It  may,  however,  be  more  convenient  to  write  F.q,  30  as  a  function  of 
('Pi  >  %)  rather  than  (xj  ,  )^),  Equation  30  represents  a  surface  over 
the  (t, ,  Af , )  plane,  the  height  of  which  measures  the  density  of 
arrivals  as  a  function  of  time  and  doppler  frequency.  If  Eq.  30  be 
written 


’d«N, } 

rdE(N,) 

dxx 

f  * 

\t(|  )  ) 

(31) 


then  if  Eq.  31  be  evaluated  at  Af*  =  Af(f  (a  constant),  there  results 
the  time  rate  of  arrivals  per  unit  doppler  frequency  at  Af».  Con¬ 
versely,  if  Eq.  31  be  evaluated  at  t*  =  t*  (a  constant),  there  results 
the  spectral  density  of  arrivals  per  unit  time  at  t^.  The  surface 
represented  by  Eq.  31  is  shown  for  typical  values  of  bottom- bounce 


25 


geometry  (Eq.  8)  in  Fig.  8.  In  this  figure,  in  (tn ,  Afu)- space  not 
enclosed  by  the  zero  contour  arrivals  are  impossible;  the  liaes  of  in¬ 
finite  arrival  density  are  a  set  of  measure  zeix)  in  (t* ,  Af* )-space, 
consequently  total  arrivals  remain  bounded  as  previous  analysis  has 
shoivn.  The  relationship  of  Figs.  4  and  7  to  Fig.  8  is  this:  Ihe 
(t,,  Af,  )-space  of  Fig.  8  is  the  transformation  frcm  (cDi ,  )- space 

of  Fig.  4  by  means  of  the  Jacobian  |j(cpi,  %)/(t, ,  Af,)!;  Fig.  7A  rep¬ 
resents  the  time  rate  of  change  of  the  cross-sectional  area  of  a  slice 
through  the  surface  of  Fig,  8  for  a  constant-time  section  and  Fig,  7B 
is  the  corresponding  change  of  cross-sectional  area  in  the  orthogonal 
direction.  Shown  also  in  Fig.  8  in  Sec.  A-A  is  a  typical  time  race 
of  arrivals  plot,  and  in  Sec.  B-B  a  typical  doppler  spectrum.  In  this 
section,  finite  beams  for  transmitting  and  receiving  may  cause  inter- 
mi  ttencies,  i.e.,  because  of  beam  limits,  certain  t,  -  Af|,  pairs  are 
prohibited. 

C.  LONG  PULSES 

If  the  pulse  transmitted  to  R  and  reflected  is  not  very  narrow, 
then  conceivably  arrivals  at  S  from  discrete  points  of  x  will  overlap 
in  time  upon  return  tc  S.  Suppose  the  pulse  transmitted  by  S  has  dura 
tion  AT.  An  arrival  at  time  t«  may  then  be  a  leading  edge  of  a  pulse 
appropriate  to  tn ,  or  to  some  portion  of  a  pulse  the  leading  edge  of 
which  occurred  prior  to  .  Thus  at  time  t*,  the  rate  of  arrival  of 
long  pulses  is 


'dECffe  )■ 

d 

r' 

‘dEClt )' 

at  '  ^  • 

A, -AT  t,i„ 

dt» 


(32) 
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whero  the  intec/rand  comes  from  Eq.  28.  The  expected  number  of  discrete 
arriv'als  at  S  at  t,  is  just  the  integral  of  Eq.  32,  i.e., 


Where  ^  ^  mean  value  of  ^  j  during  t,  -  At  to 

shows  an  extension  of  Fig.  7A  for  a  long  pulse  with 
At/Czs/tt)  =  0.004  which  for  values  of  Eq.  8  corresponds  to  At  -  10.6 

msec.  In  Fig.  9  the  maximum  of  —  0.22  jczs/n)* 


thus  for  the  typical  values  of  Eq.  8 


expected  number  of  discrete,  simultaneous  arrivals : 


i.e.,  at  t(,(zs/c)  ss  11.702,  there  are  arrivals  at  S  reflected  from  R 
along  an  expected  number  of  3.7  discrete  paths.  Referring  to  Eq.  33, 
if  it  is  desired  that  no  more  than  arrivals  are  to  be 

expected  ar  given  instant  ,  then 
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Signal  round-trtp  htintit  Hma,  t|^/(Z^/C) 


FIGURE  9  Long  Pulse  Discrete  Arrival  Oansity  vt.  Transit  Time 


Using  Fig.  7A  and  the  non-dimensional  form  of  Eq.  34  gives 


2.64  msec  pulse  length.  This  latter  calculation  gives  some  measure  to 
the  meaning  of  "narrow  pulse.”  It  appears  reasonable  to  call  a  narrow 
pulse  one  such  that  as  a  function  of  the  geometry  of  concern,  i.e., 


Fig.  4,  it  is  expected  that  no  more  than  one  round-trip  arrival  will 
occur  at  a  given  instant.  Some  intuition  to  the  relation  between 
pulse  width,  beamwidth  and  expected  number  of  discrete  arrivals  may  be 
gained  by  over-plotting  on  Fig.  5  time  bands  corresponding  to  pulse 
widths  of  AT ;  the  integral  over  the  angular  intervals  corresponding  to 
the  time  bands  yields  also  [E(14)]^^.  This  is  shown  in  Fig.  10  where 
transmitting  and  receiving  beamvddths  of  10  degrees  centered  at  a  de¬ 
pression  angle  of  20  degrees  are  Included  to  show  the  truncating  effect 
on  expected  arrivals  due  to  restricted  beams.  Apparently  from  Fig.  10, 
after  a  sufficient  period  of  time  that  an  arbitrarily  large  fraction  of 
the  pulse  width  has  passed  from  the  (cpi  ,  c|^)-space,  a  succeeding  pulse 
may  be  transmitted  with  an  expected  degree  of  ambiguity  resulting. 

Thus  a  basis  for  pulse  repetition  frequency  is  established  for  bottom- 
bounce  mode  of  propagation  of  sound. 

The  spectral  distribution  of  a  long  pulse  is,  with  reference  to 
Fig,  8,  the  integral  over  the  surface  between  (t(,  -  at)  and  t|,  and  an 


Receive  beam  ray  depression  angle,  (deg) 
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Transmit  beom  ray  depraiiion  angle,  (deg) 


FIGURE  10  Long  Pulse  Display  in  Transmit-- and  Receive— Beam  Angle  Space 


infirdtesimally  wide  sti’ip  dAf*  .  Thus,  using  Eq.  31,  the  spectral 
distribution  of  arrivals  at  t*  for  pulse  of  width  At  is 


/dECNj  )\/dE(N,  )' 

\l  T  (Xi  ,Xe  )  1 

^  )  1 

(ti,At.)| 

dt,  (35) 


Section  B-B  of  Fig.  8  provides  a  basis  for  long  pulse  return  spectra; 
the  integration  of  Sec.  B-B  over  At  will  have  a  smoothing  effect  on  the 
section.  The  resulting  spectrum  for  pulse  length  of  At  5^  0  will  prc- 
gress  from  that  shown  in  the  section,  to  one  of  roughly  Gaussian  char¬ 
acter,  then  eventually  to  a  more  or  less  flat  spectrum  for  very  long 
pulses. 
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III.  PEKTURBATIONS 


The  analysis  in  the  foregoing  admits  kinematics  and  bottom  slopes  but 
keeps  the  geometry  of  Fig.  4  static  and  assumes  a  planar  bottom  for 
computing  9|(,  and  t  and  Af.  The  error  in  this  needs  to  be  considered 
as  a  qualification  to  the  analysis.  Now  from  Eq.  2,  0„  may  be  written 

as 

cp  -  arccot  ^  cot  cp^j  (36) 

The  uncertainty  in  x,  arises  in  the  motion  of  S  and  R  during  At,  i.e., 
during  At,  x*  changes  by  an  amount 

Ax,  =  (V|  +  V,)  AT 

If  Oq  «  1,  Vs  -  V,  and  ss  z,  then 

^x,^  -  Ax,  2:  2z,  cot  %  -  2v8At 

=  2z,  cot  cp.  1^1  -  (j^)  f.] 

where  the  second  element  in  the  brackets  is  an  error  term.  For  the 
typical  values  used  previously 


(37) 


(38) 
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0.000015 


(^)(-ik)  ‘““P-  =  (o.OOl)  (0.004)  (0.364) 


is  o(lO~®) . 


i.e.,  the  fractional  error  in  x*  due  to  Vj  and  v,  during  At  is 
Another  unrelated  uncertainty  arises  if  x»  is  estimated  from  a  knowl¬ 
edge  of  cpa  and  <Zs>  ;  Xn  may  then  be  uncertain  by  approximately 
^aj/<Zs>^  cot  cfo )  i.e.,  about  400  ft  for  conditions  of  Eq.  8  (Eq.  40) 
The  uncertainty  in  Zs  and  z,  arises  from  the  fact  that  a  surface 
having  a  slope  standard  deviation  a  has  also  an  elevation  standard 
deviation  a* ,  so  that  the  distances  of  S  and  R  above  the  point  of  re¬ 


flection  are  uncertain  in  the  order  of  a.  .  In 


present  problem  we 


have  assumed  values  of  Og  and  of  [-Bg(0)/Bg(0) j  and  we  nay  approximate 
0,  from  these.*  Now,  approximately 


‘Jq  S'  o./r. 


(39) 


where  r,  is  the  correlation  length  of  bottom  elevation.  If  r,  =  OCr^) 
o|[-  Bj(0)/3Q(0)j^  I  as  is  likely,  then 


=  o{o,  [-  b;(o)/b,(o)]-*} 


(40) 


In  what  follows  some  intuitively  based  approximationo  are  used.  In 
parallel  with  0  =  arctan  dz/dx  ~  Az/Ax  if  Az  «  Ax,  the  assumption 
is  made  that  Og  ^  Of/Vi.  Also,  if  for  exa.mple  6g(r)  is  of  the  form 

Bg(r)  =  exp  j^-(r/rQ)^j,  then  |^-B"(0)/Bg(0)j"^  • 
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For  the  values  of  Eq,  8  used  previously 


and  from  Eq.  25,  the  magnitude  of  uncertainty  in  t*  is  At,/(zs/c)  ^ 

2 /<zs>)  CSC  cpo  which  has  a  typical  value  of  0.043.  A  com¬ 
parison  of  the  value  At^/zs/c  =  0.043  with  the  time  interval  of  Fig. 

7A  shows  that  with  80  percent  of  the  arrivals  occurring  within  0.045, 
importantly  large  time  dislocations  of  expected  arrivals  occur  due  to 
bottom  irregularities.  Thus,  for  the  situation  in  which  correlation 
length  for  bottom  elevation  is  very  small  compared  with  the  ensonified 
interval,  arrivals  will  occur  more  or  less  randomly  in  time  and  the 
standard  deviation  of  the  time  uncertainty  tends  to  be  a  2  ^(a,/c) 
CSC  cPo .  Hence  the  situation  described  in  Figs,  2,  3,  7A,  8,  9,  and  10 
and  the  pulse  repetition  frequency  criterion  are  ideals  to  be  hoped 
for  but  not  to  be  expected,  A  more  realistic  estimate  of  arrival  rate 
(as  in  Fig,  7A)  of  narrow  pulses  is  described  by 


’dE(lt)‘ 

1 

(41) 


in  which  is  the  standard  deviation  of  round-trip  travel  time 

2^  (cJ*/c)  CSC  cpo^  and  a  random  normal  distribution  of  z  hence  t* 

is  assumed.  The  operation  of  Eq,  41  tends  of  course  to  diminish  the 

It  is  assumed  that  at*  is  proportional  to  the  square  root  of  the  sums 
of  variances  of  z  for  each  leg  of  the  round-trip  path, 
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maximum  rate  of  arrivals  and  the  number  of  simultaneous  arrivals  for 
vd.de  pulses  and  to  distribute  arrivals  over  a  greater  interval  of 
time.  Figure  7A  contains  a  dashed  estimate  of  the  effect  of  Eq.  41 
upon  the  a*  =  0  case  calculated  there. 

The  effect  of  bottom  irregularities  upon  doppler  frequency  is 
much  less  pronounced.  Equation  26  shows  Af*  to  be  uncertain  to  the 
degree  that  9„  is  uncercain;  this  has  been  shown  to  be  very  small. 
Thus  the  doppler  frequency  distributions  of  Figs.  2,  3,  7B,  and  8  are 
satisfactorily  descriptive,  but  Fig.  8  vdll  be  smoothed  as  mentioned 
above  by  time-of -arrival  uncertainty. 
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IV.  CONCLUSIONS 


This  paper  contains  an  analysis  for  determining  the  time  of  arrival 
and  doppler  frequency  characteristics  of  a  pulse  transmitted  to  a 
target  by  reflection  from  a  not-too-rough  surface  such  as  the  sea 
bottom,  and  gives  numerical  values  for  a  typical  sonar  applicction. 
The  analysis  may  be  important  in  signal  processing  associated  with 
one-vi/ay  transmission  (communication  and  interference)  anc|  round-trip 
cransmission  (ranging)  and  in  the  interpretation  of  data  from  marine 
geophysical  surveys. 

The  analysis  shows  that  given  sonar  beam  and  source-targe  .  geom¬ 
etries,  the  expected  number  of  one-way  or  round-trip  paths  between 
the  source  and  target  is  determined  primarily  by  the  variance  and 
correlation  length  of  surface  slope,  but  that  the  incremental  time 
of  arrival  of  the  propagation  aloig  these  paths  is  determined  pri¬ 
marily  by  the  variance yof  surface  elevation  from  which  slope  variance 
do  .'Ives.  Contrasted  to  dependence  of  time  of  arrival  upon  surface 
elevation  variance,  doppler  frequency  is  dependent  primarily  upon 
surface  slope  variance.  Thus,  given  the  relationships  among  vari¬ 
ances  and  correlation  lengths  for  surface  elevation  and  some  of  its 
derivatives  (p.  33),  the  expected  number  and  doppler  frequency  char¬ 
acteristics  are  determined  primarily  by  the  slope  statistics  of  the 
surface  and  time -of -arrival  characteristics  by  surface  elevation 
statistics . 
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Furthermore,  the  analysis  shows  a  basis  for  maximizing  or  other¬ 
wise  modifying  the  expected  number  of  arrivals  for  a  given  beamwidth 
through  so^ar  beam  positioning  or  through  tailoring  of  pulse  length 
or  pulse  repetition  frequency.  A  criterion  for  pulse  length  and 
pulse  repetition  frequency  for  bottom-bounce  propagation  emerges 
from  the  analysis . 

Lastly,  the  analysis  need  not  be  restricted  to  the  sonar  problem 
described  here,  but  may  be  applicable  to  ionospheric  and  other  mete¬ 
orological  investigations  in  which  reflecting  or  refracting  centers 
may  be  in  motion  and  source  and  target  may  or  may  not  be  moving. 
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APPENDIX  I 


DERIVATION  OF  RELATION  FOR  EXPECTED  NUMBER  OF  MULTIPATHS 


Consider  a  counting  furiCtion  F4  such  tuat 


fi (e«,  e)  ^  i  when  I e„  -  e|  <  d/2 
=  0  otherwise 


I-l 


where  9  =  3(>;)  and  QmCx)  are  functions  of  x.  Then,  in  the  region  of 
X  of  interest,  i,e.,  Ax  the  numier  of  one-way  crossings  Nj  will  be 
(Ref.  4) 


11m  /  I  e'  i  Fa  (6m,  9)  Jx 

d-G  •'Ax 


1-2 


Now  i'»upi>0S6  that  9^(x)  is  a  random  process  with  W(0m  ,  9^:  x)  as  joint 
probability  density  function  with  0m*  Then  the  expected  number  of 
crossings  of  (9m,  0),  i.e.,  the  number  of  one-way  multipaths  over  a 
one-dimensionally  rough  surface  is  given  as 


E(Nx  ) 


;  e'  I  WOm, 


9  ;  x)  d9  dx 


1-3 
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Now  if  0  and  9^  have  normal  distribution,  then 


°  S.raw-ci.  f  f  l«'l  ®=<P 

9  9  *'4x  *'-» 


dx  1-4 


Integrating  first  with  respect  to  6  *  yie.  .ds 


E(Nj  )  = 


=  ^  (v)  Xx  [-  ^ 


dx 


1-5 


But  if  Bg(r)  is  the  correlation  function  of  0  and  if  Bg/(r)  is  the 
correlation  function  of  0^  then  it  is  true  Bg(0)  =  and  Bg/(0)  = 
-Bg(0)  =  a|,,  so 


In  a  similar  v'ay,  tl'ie  expactad  numbtr  of  one-way  multipaths  for 
two-dimensional  roughness  is  (Ref.  2) 


in  which  =  i>K(x,y)  and  for  round-trip  multipaths  for  two-dimensional 
roughness 
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1-8 


dxidxe  dyidyfe 


and  this  could  of  course  be  elaborated  to  round  trips  over  two- 
dimensionally  rough  surfaces  bounding  both  surface  and  bottom,  in 
which  case  an  eight-fold  integral  would  result  (Ref.  2). 

What  has  gone  before  assumes  9  is  randomly  distributed  and  z  is 
precisely  known,  which  may  not  be  true.  If  9  and  z  are  simultaneously 
uncertain  than  the  counting  function  becomes  for  z  =  z  -  e 

F.(9.  z)  =  ^  when  I  cot  (»  -  29)  -  ^  |  §  1-9 

=  0  otherwise 


In  this  event  E(Ni )  becomes  (for  z'  ^  6) 

I  cot  [cp  -  20]  - 
•  W  [z,  9,  9':  x(cp)]  dz  d0  d0'  dx 

integrated  over  the  whole  domain  of  z,  9,  0^  and  x.  The  complexity 
of  this  for  the  purpose  of  evaluatioi.  makes  tha  approximate  approach 
of  Eq.  41  appealing. 
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APPENDIX  II 


DOPPLER  FREQUENCY  RELATIONSHIPS 

Let  Fig.  II-l  be  used  to  represent  doppler  frequency  for  both  one-way 
and  round-trip  paths.  For  the  one-way  path  if  fo  is  the  at-rest  radi- 
ation  frequency  of  S  then  the  frequency  fj^  received  at  R  is 


^  =  1  +  (v,/c)  cos  (cpi  -  29i ) 

fo  1  -  (Vs/C)  cos  iPi 


II-l 


If  Vj/c  «  1  and  Vti/c  «  1,  as  is  usual,  then 


»  (V|/c)  cos  «P|  +  (V||/c)  cos  («Pi  -  26, )  II-2 

In  the  round-trip  case,  the  velocity  v*  may  be  resolved  into  a 
component  bisecting  the  arriving  and  returning  rays,  in  which  case, 
using  the  approximation  just  developed 


(vs/c)  cos  cpi  +  (Vs/o)  cos  % 


+  2(v,/c)  cos|^  [(Cft8-2ea)  +  (cpi -201  )j|cos|-5  [(Cfij-208  ) 
See,  for  example,  any  elementary  physics  text. 


II- 3 

('Pi  “201 
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Applying  suitable  trigonometric  identities  to  this  gives 

(Af(,/fo)  =  (vs/c)  cos  cpi  +  (vn/c)  cos  (cpi-29i) 

+  (vs/c)  cos  cfig  +  (vn/c)  cos  (c(:fe-202  ) 


II-4 


that  is,  for  the  round-trip  case  xt  is  as  if  a  doppler  frequency 


=  cos  cpi  +  (vn/vs)  cos  (cpi-29i),  i  =  1,2  II- 5 


may  be  associated  with  each  leg  of  the  round-trip  path. 
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13  ABSTRACT 


This  paper  contains  an  analysis  for  determining  the  time  of 
arrival  and  doppler  frequency  characteristics  of  a  pulse  transmitted 
to  a  target  by  reflection  from  a  not-too-rough  surface  such  as  the 
sea  bottom,  and  gives  numerical  values  for  a  typical  sonar  applica¬ 
tion.  The  analysis  may  be  Important  In  sonar  signal  processing 
associated  with  one-way  or  round-trip  transmission  and  In  the  Inter¬ 
pretation  cf  data  from  marine  geophysical  surveys;  It  may  In  addition 
have  application  ';o  propagation  of  electro-mechanical  v^aves  by 
ionospheric  refraction. 
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